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The anyon exciton model is generalized to the case of a neutral exciton consisting of a valence 
hole and an arbitrary number N of fractionally-charged quasielectrons (anyons). A complete set of 
exciton basis functions is obtained and these functions are classified using a result from the theory of 
partitions. Expressions are derived for the inter-particle interaction matrix elements of a six-particle 
system (N = 5), which describes an exciton against the background of an incompressible quantum 
liquid with filling factor v = 2/5. Several exact results are obtained in a boson approximation, 
including the binding energy of a (N + l)-particle exciton with zero in-plane momentum and zero 
internal angular momentum. 
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I. INTRODUCTION 

The anyon exciton model (AEM), which considers a neutral exciton made up of a valence hole and several 
fractionally-charged quasielectrons with fractional statistics (anyons), was proposed a decade ago by Rashba and 
Portnoii It was later developed^ to model incompressible quantum liquids (IQLs) with filling factors v = 1/3 and 
2/3, and is valid at large separation between the hole and the two-dimensional electron gas (2DEG), so that the 
Coulomb field of the hole does not destroy the IQL. The AEM was motivated by experiments on intrinsic photolumi- 
nescence (PL) in the fractional quantum Hall regim e) 3 : 4 : 5 that showed a pronounced double-peak structure in emission 
spectra, and it has provided major insights into the role of electron-hole separation in determining the optical spectra. 
The model yields multiple-branch energy spectra and gives a full classification of states for a four-particle system. 
It also predicts an increase in ground-state angular momentum and flattening of the corresponding dispersion curve 
with increasing layer separation. The AEM has also shown good agreement for intermediate electron-hole separation 
with numerical finite-size calculations by Apalkov et al&l 

At the time when the model was proposed, most experiments on intrinsic PL were carried out at constant carrier 
density, with the filling factor changed by varying the magnetic field. However, the AEM is valid for fixed fractional 
filling factors only. More recent experimental techniques enable the 2DEG density in the PL experiments to be 
changed using applied gate voltag ea 8 . 9 . 10 : 11 or by varying the light intensity^ Thus, it is now possible to manipulate 
the effective electron- hole separation (in units of magnetic length) while keeping the filling factor constant. So far, 
charged excitons have been observed for small electron-hole separation, but for large separations it is predicted that 
neutral anyon excitons should be observed. 

These recent developments have rekindled our interest in the theoretical treatment of spatially-separated electron- 
hole systems. In this paper we address some important mathematical features of the AEM. In Sec. [H] we derive a 
general few-particle wavefunction for a neutral exciton consisting of a valence hole and N anyons. In Sec. IIIII we 
move to consider interaction matrix elements for a six-particle anyon exciton, which should be useful for analysis 
of intrinsic PL experiments at v = 2/5. (Note that there has been some debate as to whether the ground state at 
v = 1/5 is actually an IQL - see, e.g., Ref . 113! and references therein.) In the last Section, exact results for the general 
(JV + l)-particle case are obtained and discussed. It is shown that the exciton remains bound even at large separations 
between the 2DEG and the hole. 



II. EXCITON BASIS FUNCTIONS 
A. Preliminaries 

We consider an exciton consisting of a valence hole with charge +e and N anyons with charge —e/N and statistical 
factor a. The hole and anyons reside in two different layers, separated by a distance of h magnetic lengths, and are 
subject to a magnetic field H = Hz perpendicular to their planes of confinement. Unless otherwise stated, we shall 
assume that the hole and the quasielectrons are in their corresponding lowest Landau levels. 

To simplify the description of the exciton, the hole and anyons may be considered as moving in the same plane 
with coordinates and Tj , respectively, and the the layer separation h can be introduced later when considering the 
anyon-hole interaction. 
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An exciton consisting of a hole and A anyons, all in the lowest Landau level, will have a total of N + 1 degrees of 
freedom. As the exciton is neutral, we can assign it an in-plane momentum k, which absorbs two of these degrees 
of freedom. For A ^ 2 the exciton will have A — 1 internal degrees of freedom, which results in internal quantum 
numbers and a multiple-branch energy spectrum. 



B. Derivation of basis functions 



At this stage we consider the particles as non-interacting and introduce interactions later as necessary. We can 
therefore write the (A + l)-particle Hamiltonian as 

where qt = +e and q a = — e/A are the hole and anyon charges, respectively. 
Choosing the symmetric gauge 

A = J[Hxr], (2) 
and scaling all distances with the magnetic length Ijj = (ch/eH) 1 ' 2 , we obtain 

1 /1_ A 2 A 1 (I- 1 



J = l 

where, as usual, e, H, c, and the dielectric constant are assumed equal to unity. 
We now introduce the following new coordinates (see Fig. ^ : 

D 1 ( 1 f ^ 1 V* 1 ^ 

together with the complex coordinates Q = £, X j + i£ y j. Note also the following constraint on these coordinates: 

N N 

E«i = ECf = o. (5) 

The fact that the variables £j are not independent means that derivatives with respect to these variables are 
not defined. However, it is possible to introduce the derivatives V^., which treat the variables ^ as if they were 
independent^ by working in a space of higher dimensionality and then considering a submanifold in this space 
defined by constraint ©. 

Using the chain rule, we may now rewrite the derivatives of a function / with respect to Vj as 

N 



i w 

= V €j /(...,^,...)--]T V 5 /«.... (6) 



A 
z=i 



The old coordinates can be expressed in terms of the new ones as follows: 



r h =R+|, r i =R-|+^. ! (7) 



and the corresponding derivatives as 



V, l = iv R + V p , (8) 

11 1 W 

V ^2A Vr -A V ' + V ^- A ^ V ^ ( 9 ) 
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z-axis 




FIG. 1: Two-dimensional coordinate system for hole and several quasielectrons (anyons). The hole is considered as being in 
the same plane as the anyons, with layer separation introduced when considering the anyon-hole interaction. 'C indicates the 
centre of negative charge. 



We now write the Hamiltonian J3J in terms of the new variables. Using Eqs. we obtain 

#o = ' ' ■ , , + , 

where 



(10) 



1 



2m/, I 2i 



1 



2Nm a \2i 



1 



1 



R 



i x (R+ ^ 



z x R 



and 



1 N fl 1 N 1 

^ = ^£(^.--|> ! + -[zx^ 



(11) 



(12) 



Note that Eq. (U2J) was obtained by applying constraint (J3J). 

The first part H e .xc is similar to that for a standard diamagnetic excitonji£*i& with the electron mass replaced with 
N anyon masses. The eigenfunctions of this operator can be written straightforwardly as 



*exc(R, P) = exp {ik ■ R + ii ■ [R x p] /2} $(p), 



where $(p) satisfies the equation 



1 U v [z x (p - d)] 



2nih 



1 



2Nm a [i 



1. 



[z x (p - d)] 



(13) 



(14) 



Here, d = k x z is the exciton dipole moment. Note that Eq. I|13l) does not contain the mass-dependent phase factor 
that appears in the wavefunction in Ref. This is due to our choice of the center of charge of the exciton for our 
coordinate R, rather than the center of mass chosen in Refs. ll5lH^ . Thus, the quantum number k entering Eq. i|13|) 
represents the momentum of the geometrical center of the exciton. 
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The eigenfunctions are then 

$ nm = \p- d|H £ H {(p - d) 2 /2} exp {zm0 - (p - d) 2 /4} , (15) 

where 4> is the azimuthal angle of the vector (p — d) and L™(z) are the associated Laguerre polynomials. The 
corresponding eigenvalues are 

if \m\-m + l\ if \m\+m + l\ , . 

E exc = — U+^— + — - U+J-L- . (16) 



m/, \ 2 ) Nm, 

Note that in Eqs. (|15fl and H16|l we no longer restricted ourselves to the lowest Landau level for the hole and anyons. 
For the ground-state (n,m = 0), Eq. l(T3|) reduces to 

*exc(R, p) = exp {ik • R + li • [R x p] /2 - (p - d) 2 /4} , (17) 

which depends on the quantum number k alone. 

We now move to consider the term in Eq. 1|12[) . which may be divided into two parts: 

^-^('-ffJE (7^ + ^ [»*<.]) - M> 

3=1 

and 



2« 



iSBSg { " ^ + IB I* x « J ' ( V£ . - 1 VS ') + ^ P + "* ,M) 



The eigenfunctions of the ground state of can then be written in terms of complex coordinates as 

^=Il(0) ft e x P{-IOI 2 /4iV}, (20) 



3 

with the corresponding eigenvalues 

E* = — ( 1 - — V (21) 

If we now apply constraint JSJl, we find that the function is also an eigenfunction of i?2£, with a corresponding 
eigenvalue of zero, i.e. 

# 2 ^« ^ = 0- (22) 

A calculation of the total energy of the ground state gives 

E = E exc + E i = ^— + ^— = Uw h + Nuj a ), (23) 

which is indeed the ground-state energy of the original (N + l)-particle Hamiltonian in Eq. 
The most general form for the ground-state eigenfunction of Hq is 

* (R, p, {&}) = exp {ik ■ R + ii ■ [R x p] /2 - (p - d) 2 /4} F{. . . , Q, . . •) [] ex P {~ICp| 2 /4A^} , (24) 



where F must only be a function of the complex conjugates Q. We choose F as follows to satisfy the interchange 
rules for anyons: 

F(...,C 1 ,...) = Fl(...,C 1 ,.-.)II(0-C/) q , (25) 

3<l 
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where Pl is a symmetric polynomial of degree L in the variables Note that for k — the problem has rotational 
symmetry about the z-axis, and the degree of the symmetric polynomial L is related to the exciton angular momentum 
[L z = —L — N(N — l)a/2]. 
We are now in a position to express the anyon exciton basis functions in the final form 

* (R, p. {Q}) = exp {zk • R + iz ■ [R x p] /2 - (p - d) 2 /4} P L (. . . , Q, . . .) J](C, - Q) a J] exp {-| Cp | 2 /47V} . (26) 

j<i p 

C. Symmetric polynomials 

We now consider the precise structure of the symmetric polynomials Pl which appear in Eq. I|2()l) . To determine the 
symmetric basis polynomials of order L we apply the fundamental theorem of symmetric polynomials^ which states 
that any symmetric polynomial in N variables, P[x\, . . . , a; at), can be uniquely expressed in terms of the elementary 
symmetric polynomials 

Q~l = S/Xj, 0~2 = / J XjXj, er 3 = ^ XiXjXk, ■■■ , cr N — XiX 2 ■ ■ ■ x N - (27) 

i i<j i<j<k 

It is now apparent that all linearly-independent symmetric basis polynomials of a particular degree may be enumerated 
by considering the possible products of the elementary symmetric polynomials (|27|l , as the total degree L is the sum 
of the degrees of the constituent polynomials. For example, the possible symmetric polynomials of degree four are 
constructed from af, cr\a2, (T10-3, erf, and 04. 

There are two further problems. The first is to determine the number of possible products of elementary symmetric 
polynomials for a total degree L. The second is to determine the structure of these products. 

To calculate the number of linearly-independent symmetric basis polynomials of degree L which may be constructed 
from the elementary symmetric polynomials ax, ... , o~n we use a result from the theory of partitions^ 

The number of ways of partitioning a number L into parts of size 1, 2, . . . , P is given by the coefficient of x L in the 
expansion of 

■ • • • 7— ^ = n — ^ ( 2g ) 

k=l 

where each term in the product is known as a generating function. The number of ways of partitioning increases rapidly 
with L and does not follow any pattern. Furthermore, the different products of elementary symmetric polynomials 
for a particular value of L must be determined by hand. An approximate formula for calculating the number of ways 
of partitioning a number does exist, the so-called Hardy-Ramanujan formula, M but it is not applicable to the present 
case as we have no polynomial of order one due to constraint JSJ) . This constraint significantly reduces the number of 
possible symmetric polynomials by removing the first factor in the product (|28J) . For example, for N — 3 the number 
of polynomials of degree L is equal to the integer part of L/6 + 1 for even L and the integer part of (L — 3)/6 + 1 for 
odd L, which corresponds to the result obtained in Ref. |2j. 

D. Discussion 

The key difference between the current general formulation and that outlined in Ref. |2j for the four-particle case 
(N = 3) is the replacement of anyon difference coordinates Zji by the new coordinates Q. The principal advantage 
of the difference coordinates was that they simplified the calculation of inter-anyon repulsion matrix elements, as 
well as the form of the statistical factor in the exciton wavefunction. However, they had the disadvantage that the 
classification of symmetric polynomials was more difficult, as it was necessary to introduce a Vandermonde determinant 
for odd-L polynomials. For example, in the simplest case L = 3, even though the number of anyon coordinates is the 
same in both formulations, we have 03 = (212 — ^23)(z23 — 231X231 — ^12) in terms of difference coordinates, whereas 
we have a simple product cr 3 = C1C2C3 m terms of the new coordinates. For a number of anyons greater than three this 
disadvantage is crucial, as the classification of polynomials in terms of Zji becomes too cumbersome and the number 
of constraints on these coordinates is greater than one. 
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TABLE I: Possible ways of constructing a symmetric polynomial Pl from the elementary symmetric polynomials 02, 03, 04, 
and <r 5 . 



Order, L 


No. of polynomials 


Structure 





1 


1 


1 





- 


2 


1 


(72 


3 


1 


CT3 


4 


2 


2 

(72, (74 


5 


2 


(72(73, (75 


6 


3 


3 2 
(72, (72(74, (73 


7 


3 


2 

<72(73, (72(75, 0"3(74 


8 


5 


4 2 2 2 
(72, (7 2 (74, (72(7 3 , (73(75, <T 4 


9 


5 


3 2 3 
(7 2 (73, CT 2 0"5, (72(73(74, (73, (74(75 


10 


7 


Crf, (72(74, Crfaf, (72(73(75, (7204, (73(74, O5 


11 


7 


4 3 2 3 
CT 20"3, 0"2°"5, (7 2 (73(74, (72(73, (72(74(75, 

0~f(75, a3Cr| 


12 


10 


Crf, Cr|(74, (72(7|, (72(73(75, Ct|ct|, CT2C7|(74, 

0"2Cr|, CT 3 , (73(74(75, crf 



III. SIX-PARTICLE ANYON EXCITON: BOSON APPROXIMATION 



A. Formulation 



We introduce an anyon exciton consisting of a hole and five anyons with charge — e/5. We also make a boson 
approximation so that the statistical factor a = 0. Our justification for this step is as follows. It was shown in Ref. 
that for large values of h (which is required for the AEM to be valid) the statistical factor a becomes unimportant, 
and the results for a — were very similar to those for a — ±1/3. We would therefore also expect this to be true 
for a — ±1/N, as this is even closer to zero. From now on we consider a boson approximation (cv = 0) for anyon 
statistics. 

The explicit form of the exciton basis functions is 

5 

*L,M,k (R, P, {Ci}) = exp {ik ■ R + ii ■ [R x p] /2 - (p - d) 2 /4} P L , M fa, >&) J[ exp { — 10| 2 /20} , (29) 

i=i 

where M enumerates different linearly-independent symmetric polynomials of degree L. 

To construct the symmetric polynomials Pl,m we need to consider only the elementary symmetric polynomials <J2, 
(73, (74, and 175. Note that o\ — because of constraint From the above, we find that the number of possible ways 
of constructing a polynomial of degree L is therefore the coefficient of x L in the expansion 
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n jh- < 30 » 



1 — x 2 1 — x 3 1 — a; 4 1 — X' 

The possible ways of constructing the first twelve polynomials are shown explicitly in Table HJ It can be seen that the 
number and complexity of the polynomials increases rapidly with the degree L. 

It is evident from Eq. (|29|l that basis functions with different values of k are orthogonal. We now expand a general 
exciton wavefunction for given k in terms of the complete set of basis functions l|29|) : 

* = !>»•*»•. (31) 



We shall show in Sec. IIIIBl that functions with different L are orthogonal. However, basis functions with the same 
value of L but different M are not necessarily orthogonal, so their scalar products (L, M\L, M') will be non-zero. We 
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therefore write the Schrodinger equation in matrix form as 

H X = £B X , (32) 

where B is the block-diagonal matrix of scalar products (the overlap matrix) and e is an energy eigenvalue. The size 
of each block in B depends on the number of different wavefunctions with given L. For example, in the case of a 
six-particle exciton the block corresponding to L = 12 will be of size 10 x 10 (see Table HJ. Both the overlap matrix 
B and the Hamiltonian matrix H are diagonal in k. Note that for k = the matrix H takes the same block-diagonal 
form as B (as will be shown in Sees. IIII ~C\ and UlI D|) . and as a result the problem becomes exactly soluble. 

We shall now proceed to evaluate the matrix elements in Eq. I|32f) . As H and B are diagonal in k, in what follows 
we consider only the matrix elements diagonal in k. Since all terms in the polynomial Pl,m are of the form of a 
product of the coordinates Q, we use the following functions in monomials in the matrix element calculations: 

5 

w , k (R, p, {&}) = Cexp {ik ■ R + iz ■ [R x p] /2 - (p - d) 2 /4} CC^CCC? 5 II ex P {-|0| 2 /20> , (33) 

3=1 

where {n} denotes the set of quantum numbers n\ to 725, and the constant C will be defined in Sec. IIII 51 The basis 
functions in terms of symmetric polynomials of degree L can be obtained as a linear combination of the monomial 
functions (|33f) . and ri\ + n 2 + n 3 + n 4 + n 5 = L. 

We also take into account constraint (JSJ via the following transformation: 

f di 

<J(€l+C2 + €3+€4 + € B ) = J 7^)2 ex P{* f -(£l+£2+€ 3 + £4 + £ 5 )}- (34) 

B. Overlap matrix elements 

To calculate the overlap matrix B, let us consider the scalar product of two monomial functions: 

(M I {«'}) = J dR J d pJ / ^1^2^3^4^5*{n }l k*K},k exp {if • (£ x + £2 + £3 + «4 + « 5 )> • (35) 

The integration over R and p gives a factor of 27rA, and this leaves the following: 

({n} I {n'}) = C 2 (2nA) J f[ M„.„, (f), (36) 

where 

M w (f) = J d£C n C'e- e/10+ ^. (37) 

If we now make the substitution £ = £e l ^ s , and let <j> be the angle between f and £, i.e. 

4> = - <Pf, (38) 

we obtain 

p2n poo 

M w (f) = e l(m " m '^ f / #/ dee 1+m+m 'exp{i(m-TO')^ + */Ccos0-e 2 /lO}. (39) 



We can now make use of Bessel's integral 

{■2ti 



Jo 

20 



d ^ e ±i(m-m')4>+iM cos <f> = 2m\ m - m '\ J\ m -m-\{ft), (40) 

J 

to reduce Eq. to the form 

M w (f) - 2 7 ril m - m 'le i (" l - m ')^7Ww(/), (41) 
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where 

POO 

M mm ,{f)= d^e +m+m 'e- e ' w J\ m -m'\{m- (42) 
Jo 

This can also be expressed^! in terms of a confluent hypergeometric function $(/?, 7; z) as 

2 (m+m')/2 5 (m+m')/2+l r ( max { m m '| + 1) 

.M w / = ; jr. — 1 ^" l/2 $ max{m, to'} + 1, m - m'\ + 1; -t) , 43 

\m — my. 

where t = 5/ 2 /2 and max{m, to'} is the largest of the integers to and to'. Eq. (|43[) can be simplified further if we 
apply the Kummer transformation 

$G8,7;-t) = e-*#(7-/3, 7 ;i), (44) 

and note that (7 — (3) is always a non-positive integer: 

7 — /3 = |m — to I — max{m, m'}. (45) 

This means that $(7 — /3, 7; t) reduces to a polynomial, and 7; — t) is then just a polynomial in t multiplied by 
e l . 

If we now perform the integration over cj)f we obtain the final form of the overlap matrix elements 

poo 5 

({n}\{n'}) = C 2 A(2n) 5 5 LL , <{{ " M„ ,/;/;. (46) 

where L = n\ + 712 + + 714 + 72,5. Note that nothing depends on k in this expression. It is convenient to define the 
constant C so that for L = L' = the matrix element ({n} \ {n 1 }) = (0|0) = 1. This yields 

° 2 = 125^- (4?) 

In the above formulation we have only considered matrix elements in terms of monomial functions. As mentioned 
earlier, matrix elements in terms of symmetric polynomials can be easily constructed as a linear combination of the 
monomial matrix elements (|46|) . 

C. Anyon-anyon interaction 

The anyon-anyon interaction has the form 

We shall only calculate the matrix elements for the first term V12 as the others follow by analogy. We begin by taking 
the Fourier transform of Eq. (|48Jl : 

Via (€1, - ^ / ^ exp {zq • (^ - g 2 )} . (49) 
This gives a matrix element in terms of monomials 



(M |Vl2| K}> = ±JdRj dp I / ^ / ^1^2^3^4^5 



(2tt) 2 7 27rg 

x *{„},k*{n'},k exp {if ■ (€1 + € 2 + £3 + £4 + £5)} ex P ' (€1 ~ £2)} • (5°) 
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Following the procedure in Sec. IIII 51 yields 

({n} \V 12 \{n'}) = ° 2A ^ f II fdffdq exp {i<j> + {m - n[) + #_(n 2 - n' 2 )} 

3=1 J J 

5 

x exp {i(j)f [ (n 3 - n' 3 ) + (n 4 - n'^) + (n 5 - n' 5 )}}M nin[ (f+)M n2n ,(f-)l[M njn ,(f), (51) 

i=3 

where A4 m m' are as defined in Eq. (|42|) . f± — f ± q, and (f>± are the phases of f±. 

We now seek to eliminate <pf by the change of variables <j) — (f>f — </> q and ip± — <p± — <pf. After substitution and 
integration over (/> q we obtain 



fy2A(n \4 J roo roc rJir 

<W|F 12 | {n'}> = ±-^5 LL ,l[^^ / 4f// dq # exp{^(L-L')} 

zo . =1 Jo Jo Jo 

5 

x exp{ii/>+(ni - n[) + iip-(n 2 - n' 2 )} M nini (/ + )M„ 2 „i (/_) M„ 3 . n / , (52) 

where L = rai + n 2 + 713 + ri4 + 715, and t/>± and /± can be expressed in terms of the variables of integration as 

e i *± = (f±qe- i *)/f ± , f ± = f + q 2 ± 2/gcos0. (53) 

It is evident from Eq. I|52[) that the anyon-anyon interaction matrix has the same block-diagonal structure as the 
overlap matrix. The integrand in Eq. I|52|l can be further reduced to a product of a polynomial in q, /, and e ±4 *, 
and an exponential factor exp(— 25/ 2 /2 — 5g 2 ). This integral can therefore be evaluated analytically for any {n} and 
{n 1 }. For the case of L = L' = 0, using Eq. I|47|) for C 2 , the anyon-anyon interaction matrix element reduces to 

(0|Ka|0) = (54) 

D. Anyon-hole interaction 

The anyon-hole interaction takes the form 



j=l j=l I 3 n < 



where 

r jh =£j~ P+ hi. (56) 
Considering only V\h, we take the Fourier transform: 



y^^^^-lf J^ Vah ^ exp ^ q ■ " P)> . ( 57 ) 



where 



9ir 

V ah (q) = —e-« h . (58) 

q 



The matrix elements are then 

(W \Vi h \ {n'}) = ~jdB.Jdpj J Av*h{q) J d^ 2 d^ 3 d^ 5 

x * { „} : k*{„'},k exp {if • (£ x + £ 2 + £ 3 + £ 4 + £ s )} exp {iq ■ - p)} . (59) 
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Integration over R and p gives 2ttA exp [—q 2 /2 — id ■ q). Then, following the procedure in Section fill Bl once again, 
we obtain 



<{n} \V lh \ {„'}) = y 2A ^) 2 fl ito-n'J f di f dq Vah {q)e-f 

3 = 1 J J 



5 

x exp - n[) + ifa [(L - L') - (m - ni)]}M ni „ ; (/+) ~[[M njn >.{f), (60) 

where M mm ' are as defined in Eq. f + = f + q, and + is the phase of f+. As before, L = n\ + n 2 + n 3 + n 4 + n 5 . 

We now eliminate fa by substituting 4> = fa — 4>q and ip = <p + — fa . Choosing the x-axis along d, the integration 
over </> q is then 

2-7T 

d0 q exp{i(j)q(L - L') - idqcosfiq} = 2ir(-i)^ L ~ L 1 J\ L - L '\(dq), (61) 

o 

which gives for the matrix elements 

({n} \V lh \ KB = _ g W^ H) l*-*'l fl r J\L-V\{dq) 

x dff d<f> exp{ifaL-L')+i^( ni -n[)}M nin{ (f+)]lM njn '(f), (62) 
Jo Jo - =2 

where ip and /+ can be expressed in terms of the variables of integration as 

e i + = (f + qe- i +)/f + , fl = f + q 2 + 2fq cos 0. (63) 

For d = 0, the Bessel function J\ L _ L ,\ (dq) entering the integrand of Eq. I|62fl is non-zero only if L = L'. This indicates 
that for the case of zero in-plane momentum {k — d — 0), the anyon-holc interaction matrix has the same block- 
diagonal structure as the anyon-anyon and overlap matrices. Therefore, states with different angular momentum L z 
decouple. 

Eq. (|62f> may be simplified further, and this yields the following expression for the anyon-hole interaction matrix 
elements for basis functions in terms of symmetric polynomials: 

/>OG 

{L,M\V ah \L',M') = - / dq exv{-bq 2 /2- qh)J\ L _ L ,\{dq)Q LM .L-M>{q), (64) 
Jo 

where Qlm,L'm'(c[) is a polynomial in q (the lowest order polynomial is Qoo,oo = 1). A similar expression to Eq. 1641) 
is obtained for a four-particle exciton in Ref. 0- For h = 0, Eq. (|64[) further reduces to an expression in terms of 
elementary functions. For the case of k = d = and L = L' — 0, the anyon-holc interaction matrix element is then 

<0|V^|0) = -^. (65) 

E. Exact results for k = 0, L = 

All the above results are simplified significantly for the state with k = and L = 0. For such a state, the (N + 1)- 
particle wavefunction (j26(l with a = reduces to the form 



* (R, p, {Q}) = C N exp {ii ■ [R x p] 12 - p 2 /4} exp j- £ |C P | 2 /4A j 



(66) 



where the normalization constant Cjv can be determined from the integral 



■jdnj J dp e-" 2 / 2 J d^ ■ ■ -d^exp • f - £ 2 



fl / </R / -^r / <lp < ~ 2 I'ltr- ■</£x^l><iEK.- f -4 2 /2Ar) 1 = ^(2^ = 1. (67) 
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For the case of k = and L = 0, the interaction matrix elements of a (N + l)-particle exciton are straightforward 
to evaluate. The anyon-anyon matrix element is 



N 



(68) 



x exp {£ ■ f - tllW) + iq ■ ({, - « 2 )J = 
For non-zero inter-plane separation h the anyon-hole matrix element can be reduced to 

x y d€i • • • d£ N exp |^ (i^ p ■ f - $/2N) + • (€ x - p) - gftj 

exp(-iVq 2 /2 - g/i) = - W-^ ef^erfc (h/Vw) , (69) 



where erfc(a;) is the complementary error function. Using the asymptotic expansion of erfc(x) it can be easily seen 
from Eq. H69f) that (0|V^/j|0) — > — 1/h as /i — * oo, as expected. 

Eqs. (|68|1 and l|69|l also allow us to calculate the critical inter-plane separation h c at which the k = 0, L = state 
becomes unbound, i.e. when 

where x c = h c /^/2N . So, for N — 3 the critical separation h c w 5.39 for N — 5 we find that /i c ~ 5.59 In , and 
for iV > 1 we have /i c w 1.32V2iV Z^f . Notably, these critical separations are well inside the region for which the 
AEM is applicable. It should be emphasized that the state with L = is not the ground state for the anyon exciton 
at large separation h. For example, for a four-particle exciton^ the ground states for large separation satisfy a 
superselection rule L = 3m, where m is an integer, and when h — * oo the ground state energy tends to its classical 
value, e= -(2/3)^ 2 /h. 

Despite the non-applicability of our model to a real physical situation at small layer separations (h < Ih), the model 
remains soluble for all values of h, including h = 0. Moreover, it has been shown in Ref.Qthat the ground state of the 
four-particle problem (N = 3) at h = is the state with k = and L — 0. We expect the same to be true for N ^ 5, 
since the anyon-hole attraction will always overcome the anyon-anyon repulsion at small inter-particle separations. It 
can be shown (in a similar way to that in Ref. 0) that the smallest average inter-particle separation corresponds to 
the L = case. Therefore, the case of k = 0, L — is of special interest since it predicts the anyon exciton binding 
energies in the limit h — > 0. 

We are now in a position to write down a general expression for the binding energy of a (N + l)-particlc exciton 
at zero 2DEG-hole separation (h = 0): 



E h = - 



((0\V ah \0) + (0\V aa \0)) = 



2V2N 



k- < 71 > 



Eq. 171|) has been written in this particular form to emphasize the key result that for any value of N there always 
exists at least one bound state of a neutral (N + l)-particle anyon exciton at h = 0. For N = 1, Eq. (|71|) yields the 
value yj it /2, which corresponds to that obtained in Ref. ^3 for a standard diamagnetic exciton. For a four-particle 
exciton (N = 3), E b = (l - \/2/6) which agrees with the result of Ref.0 Finally, for N = 5 we have a binding 

energy of (1 — V2/5)\/7r/10, which can also be obtained from Eqs. |QJ| and l|6"5*|) . 



IV. CONCLUSIONS 



The anyon exciton model has been generalized to the case of an arbitrary number of anyons and several important 
mathematical results have been obtained. Starting from the Hamiltonian for a non-interacting system of N anyons 
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and a valence hole in a quantizing magnetic field, we have obtained a complete set of exciton basis functions. These 
functions have been fully classified using a result from the theory of partitions. We have derived expressions for the 
overlap and interaction matrix elements for a six-particle system, which describes an exciton against the background 
of an IQL with filling factor v = 2/5. In the particular case of k = and L = 0, we have found an expression for the 
binding energy of a (N + l)-particle exciton, which agrees with known results for a standard diamagnetic exciton and 
a four-particle anyon exciton. We have also shown that the (N + l)-particle exciton remains bound for 2DEG-holc 
separations exceeding several magnetic lengths, when the anyon exciton model becomes applicable to real physical 
systems. 



* m.e.portnoi@e x"ac.uk| Also at A.F. Ioffe Physico- Technical Institute, St. Petersburg, Russia. 
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